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ABSTRACT 

We present a model for the rotational evolution of a young, solar-mass star interacting magnetically 
with an accretion disk. As in a previous paper (Paper I), the model includes changes in the star's 
mass and radius as it descends the Hayashi track, a decreasing accretion rate, and a prescription for 
the angular momentum transfer between the star and disk. Paper I concluded that, for the relatively 
strong magnetic coupling expected in real systems, additional processes are necessary to explain the 
existence of slowly rotating pre-main-sequence stars. In the present paper, we extend the stellar spin 
model to include the effect of a spin-down torque that arises from an accretion-powered stellar wind. 
For a range of magnetic field strengths, accretion rates, initial spin rates, and mass outflow rates, the 
modeled stars exhibit rotation periods within the range of 1-10 days in the age range of 1-3 Myr. This 
range coincides with the bulk of the observed rotation periods, with the slow rotators corresponding 
to stars with the lowest accretion rates, strongest magnetic fields, and/or highest stellar wind mass 
outflow rates. We also make a direct, quantitative comparison between the accretion-powered stellar 
wind scenario and the two types of disk-locking models (namely the X-wind and Ghosh & Lamb type 
models) and identify some remaining theoretical issues for understanding young star spins. 

Subject headings: Accretion, accretion disks, Stars: evolution. Stars: magnetic field. Stars: pre-main- 
sequence, Stars: rotation, Stars: winds, outflows 
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1. INTRODUCTION 

There remain a number of unexplained phenomena re- 
lated to the spin rates of young, low mass (< 2M(t,) stars 
(e g iHerbst et all [20071: Hcholz 2009; Ir win fc Bou^ 
mim . At the youngest stages, the fundamental dilemma 
is that a large fraction of stars rotate more slowly than 
expected, given that the stars are in the process of 
contracting and are often accreting material with high 
specific angular momentum (i.e., from a circumstellar 
disk). In the absence of significant spin-down torques, 
one would expect these stars to spin near breakup speed. 
However, observations reveal that approximately half of 
solar mass stars with ages of less than a few million 
years rota te at less than ^1 % of their breakup veloc- 
ities (e.g., IRebuh et al.l[200l IHerbst et aL l lMl IScholzl 
I2Q09D . Thus, for a large fraction of pre-main-sequence 
stars, there exists some mechanism that removes signifi- 
cant amounts of angula r momentum. 

In a previous paper (jMatt et al.l 120101 hereafter Pa- 
per I), we examined whether the magnetic connection 
between a young star and surrounding accretion disk 
can be solely responsible for the angular momentum 
loss. That work employed a model similar to some previ- 
ou s stud i es in the literature (iCameron fc C ampbell 1993; 
Yil[T99l [19951 ICameron et al.l 1199,1 lAmhtage fc Clarkei 



1996( ). which computes the time-evolution of the mass, 
radius, and rotation rate of an accreting star. Specifi- 
cally, Paper I presented calculations of the spin rate of 
a one-solar-mass star as it evolves from an age of 30,000 
yr to 3Myr. The calcu lations consider ed the contraction 
of the star during the iHavashil ()1961h phase and an ex- 



ponentially decreasing accretion rate. The calculations 
also considered a range of initial accretion rates, spin 
rates, and magnetic field strengths, approximately rep- 
resentative of typical observed or expected values for low 
mass pre-main-sequence stars. Rather than attempt to 
explain all phenomena related to young star spins, the 
primary goal in Paper I (and in the present paper) was 
to address the fundamental question of whether the mod- 
els could produce spin rates within the observed typical 
range of ~ 1-10 days. 

The main difference with previous works is that the 
model in Pap er I used the m agnetic torque formulation 
presented by iMatt fc Pudritz (2005b), which includes 
the effects of the magnetic coupling (diffusivity) in the 
disk and the loss of magnetic connection that occurs 
when field lines are su fficiently twisted azimuthally (e.g., 
lUzde nskv et al.l |2002[ ). When considering the strong 
magnetic coupling (i.e., low diffusivity or high magnetic 
Reynolds number) expected in real systems. Paper I con- 
cluded that the spin-down torques arising from the star- 
disk interaction alone were not sufficient to produce ro- 
tation periods longer than ~ 3 days at ages of 1-3 Myr. 
The main result of Paper I suggests that additional spin- 
down torques are necessary to explain the slow rotators. 
That model neglected any torques that may arise along 
open field regions, such as torques from stellar winds, so 
a natural next step is to consider whether or under which 
circumstances such torques may be important. 

Powerful and large-scale jets and outflows are ob- 
served to emanate from accreting young stars with 
outflow rates of the order of 0.01-0.1 times the accretion 



2 



Matt et al. 



rates fe.g- IReipurth fc Ballvl[200llCabrid[2007l) . Stellar 
winds — that is, outflows that are magnetically connected 
to the star, as opposed to the accretion disk — may be 
an import ant co mponent in these large-sca l e flows (e.g., 
Decampli' 'igSl"; 'Kwan & Tadcmaru '1988'; 'Fendt ct aL| 

Paatz & Camc nzinc 
" 200C 



Canienzind 1996, 



1995: Fendt _ 

1991 iHirose et alJ liggTtTlp^ndt fc Elstneii 
Bogovalov fc TsinganosI l2001t ISautv et aLl 120021 



20091: iRomanova et alJ 120091 : iSautv et al 
thermore, high-resolution spectroscopic 



Meliani et alJ 120061: iMatsakos et alJ |2008[ [2O09: 

[201 




201 ID . 

observations 



of these systems show some emmision line features 
that appear to be best explained by the presence of 
stellar winds, and that these signatures generally corre- 
late wit h the accretion rates (e.g., Hartman n et al 
iKwan fc Tadcmaru 



1982; 



199C; 



1988; Hartmann et al 



jEdwards ct al. 120031: iPu prcc ct al^ 
Edwards et al.. ,2006; .Kurosawa et al.i ,2006; Kwan et al 



2005 



2007 



a ei ai.. ,zuuu; .ivmuaawa ei ai.i,zuu»: ivwau ei ai 
Johns-Krull fc Herczeej 120071: Kurosawa et al 



20111) . The idea that powerful stellar winds could also 



be the primary agent for removing angular momen- 
tum from the star has been suggested or explored by 
several authors (iHa rtmann fc MacGrcgor 1982; Mcste l 
198^ iHartmann fc Stauffcr 198 9t " pTout fc PringlS 

ijHh 



1992: iPaatz fc Canienzind il996i: iHirose et al 
Ferreira et al.l 200Ct iRekowski fc Brandenburg 



Romanova et al.' 20(1 
Rekowski fc Brandcnbur; 



.M att fc Pudritz 

2008bl: IRomanova et alT 2009: Zanni fc Ferreira 200 



iMatt fc Pud rit^ 
20061: ' 




Sautv et al.l I2011D . IMatt fc Pudritz (2005a, 2008 



showed that a stellar wind will be much more effective 
at spinning down a star than the magnetic connection 
between the star and disk, as long as the mass outflow 
rate is high enough. In order to maintain a sufficiently 
high mass loss rate, ,Matt fc Pud ritz (2005.a) suggested 
that the stellar win d is somehow powered by the ac- 
cretion process, and IMatt fc Pudritz! ()2008aD suggested 
that accretion onto the star may power a stellar wind 
by the excitation of a large flux of Alfven waves along 
the open field lines. The idea of a wave-driven and 
accretion - powered ste llar wind was further explored by 
ICranmeH(|2008l . l2009h . 

In the present paper, we explore the role of stellar 
winds in the spin evolution of a young star that is also 
interacting magnetically with an accretion disk. The 
model (described in ^J2]) is identical to that of Paper I ex- 
cept that the star experiences an additional torque from 
a stellar wind. This introduces only one new parame- 
ter, the mass outflow rate in the wind, which we take 
to be a fraction of the accretion rate, as expected for an 
accretion-powered stellar wind (hereafter, APSW). The 
results in section [3] demonstrate how an APSW, acting 
during the entire accretion history, may explain the typi- 
cal range of observed spin rates. In section|4l we compare 
the present model with disk-locking models. A discussion 
of the conclusions of this work and remaining challenges 
are contained in section [5] 

2. STELLAR SPIN EVOLUTION MODEL 

In order to best compare the effects of powerful stel- 
lar winds to the torques arising in the star-disk mag- 
netic interaction alone, we employ the same assumptions, 
equations, and method as Paper I, except that we have 
included the mass loss and torque resulting from the 



wind. This section lists the relevant equations and points 
out where there are differences from Paper I (where the 
reader will also find more details). 

2.1. Mass Flow 

The mass accretion rate onto the stellar surface follows 
an exponential decay, 

Mao e-*/*% (1) 



Ma 



where Alao is the mass loss rate a.t t — 0, and ta is the 
exponential decay timescale. Following Paper I, we adopt 
ta = 10^ yr, for all models, and two different values of 
Alao so that the models span a range of accretion rates. 
The two values correspond to a "low" accretion rate with 
Mao ~ 1O~^M0 yr"-'^ and a "high" accretion rate with 

Mao = 10-'^ Mq yr-i (see Fig. 1 of Paper £])• In ah cases, 
we chose an initial stellar mass, such that the mass of the 
star at the end of the computation (at an age of 3 Myr) 
equals one solar mass (see Fig. 2 of Paper I). 

At the same time, the present work considers the ad- 
ditional effect of a stellar wind. The general picture of 
simultaneous infiow and outflow from the star is as fol- 
lows. A fraction of the stellar magnetic flux connects to 
the accretion disk, which is responsible for the magnetic 
torques between the two and for channeling an accretion 
flow onto the star. The stellar wind fiows from the polar 
regions of the star, where the magnetic field lines reach 
beyond the magnetosphere. This region of the magnetic 
field is susceptible to being opened and kept open by 
energ etic processes, such as disk and/o r stellar winds 
(e.g., ISafiCTllT99l IMatt fc Pud riti l2005bl) or differential 
rotat ion between the star and disk (e.g.. lUzdenskv et all 
I2002f) . The basic idea behind APSWs is that a fraction 
of the accretion energy that is dissipated near the sur- 
face of the star ultimately powers a stellar wind. For this 
to work, energy deposited along the accreting field lines 
should transfer laterally (likely via waves) across the ac- 
cretion flow and into the open field region. In the open 
field region, this energy then drives an enhanced stellar 
wind, where the outfiowing mass is from the star and is 
distinct from the material that is currently accreting onto 
the star. A more detailed descript i on of t hese processes 
are gi ven m (e.g.) IMatt fc Pudritd (|2005aD and fCranmeil 
(|2008f ). and an illustration of the model is given in Figure 
1 of each of those works. 

When the energy responsible for driving the wind ulti- 
mately comes from the accretion process, it stands to rea- 
son that the mass outfl ow rate Mw will be ultimately tied 
to the accr etion rate ()Matt fc Pudritzl [2005bt iCranmeJ 
120081 [2OO9I) . Thus, we parameterize the stellar wind out- 
flow rate using 



X ■ 



The evolution of the star's mass is thus given by 



dt 



(1 - X)Ma. 



(2) 



(3) 



^ Paper I employed a different nomenclature, where MaO = 
Mjyta^ , but regardless we adopt the same two values for the initial 
accretion rates as in Paper I. 



Spin Evolution: Effect of Accretion-Powered Stellar Winds 



3 



iHartmann fc Stauffeii (|1989[ ) suggested that a value of 
X ~ 0.1 is necessary for the angular momentum outflow 
rate in a stellar wind to be comparable to the angular mo- 
mentum deposite d onto the star from ac cretion of disk 
material (and see lMatt fc Pudritzl 1200 5b'). Observations 
of large-scale outflows from T Tauri stars and (presum- 
ably younger) Class I sources indicate mass outflow rates 
of ~ 0.1 times the accretion rates, with about an order 
of magnitude of sc atter and/or uncertainty in this value 
e.g., see review bv lCabrit]|2007D . It is still not clear what 
fraction of the mass outflow rates observed on large-scales 
may be due to flows that are magnetically connected to 
the star (as required in the present work), as opposed to 
flows from the accretion disk (e.g.. iFerreira et al.l 120061: 
Kurosawa et al.|[2006l : iKwan et al.|[2007HKurosawa et al.l 



20111). Recent models of accretion-powered stellar winds 
bv ICranme'^ (|2008l 120091 ) exhibit outflows from the star 
itself with x 0.01. Thus, in the present work, we con- 
sider cases with both x — 0.1 and 0.01, to explore a 
possible range of this value. 

2.2. Stellar Structure and Evolution 

The evolution of the stellar structure follows simple 
Kelvin-Helmholz contraction of a polytrope (with index 
3/2). This treatment results in an evolution of stellar 
radius that follows 



^ = 2— M.(l-x) 



28naRtTf 
3GM; 



(4) 



where i?, and are the stellar radius and mass, 
is the effective temperature of the star, and a and G 
are the Stefan-Boltzmann constant and Newton's grav- 
itational constant. As in Paper I, we adopt an initial 
stellar radius of 8Rq and constant photospheric tem- 
perature Te = 4280 K for all calculations, so that the 
evolution of the star 's struc ture resembles that from the 
model of iSiess et al.l ()2000[ ). as the star descends along 
the Hayashi track in the H-R diagram (see Fig. 3 in Pa- 
per I) . The only difference between this equation and the 
corresponding equation of Paper I is the addition of the 
factor (1 — x), which takes mass loss into account. 

2.3. Spin Evolution 

The evolution of the angular spin rate of the star fol- 
lows 

dn, / Ma,^ , , 2 dRA 

-w-i:-''*[m:^'-''^-^r:^)^ 

where fi* is the (solid body) angular rotation rate of the 
star, is the net torque on the star, and /* is the stellar 
moment of inertia, /, = fc^A/»_Rj, where k is the nor- 
malized radius of gyration (we adopt = 0.2). The only 
difference between this equation and the corresponding 
equation of Paper I is the addition of the factor (1 — x), 
which takes mass loss into account. 

It is often instructive to express the spin rate as a frac- 
tion of the breakup speed, defined as the Keplerian ve- 
locity at the star's equator. This normalized spin rate is 
defined 
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Following Paper I, we consider two cases with different 
initial spin rates. The two cases have initial fractional 
spin rates of /o = 0.3 and 0.06, representing the two 
extremes of rapid and slow initial rotation. 

2.4. Torques on the Star 

The present work considers the simultaneous effects 
of torques arising in the magnetic star-disk interaction 
and angular momentum loss from stellar winds. In both 
cases, the angular momentum gain and loss of the star is 
primarily transmitted by a magnetic field. The model as- 
sumes a rotation-axis-aligned dipole magnetic field, with 
a strength of at the surface and equator of the star. 
We consider two different values of the magnetic field 
strength, i?, = 500 G and 2000 G, in addition to a case 
with i?* = used for comparison. The model assumes 
that is constant in time for all cases. 

2.4.1. Magnetic Star-Disk Interaction Torques 

The torques arising from the magneti c interaction be- 
tween the star and accretion disk follow [Matt fc Pudritd 
(r2005b), and the governing equations are given in Pa- 
per I. It is convenient to define the corotation radius 



Rrn = 



1/3 



(6) 



which is the singular radius in the Keplerian disk that 
rotates at the same angular speed as the star. This lo- 
cation is physically meaningful in the magnetic star-disk 
interaction because it is where the differential rotation 
between the star and disk equals zero. 

The magnetic field of the star is typically strong 
enough to be able to truncate the disk at a distance of 
a few stellar radii. The truncation radius is denoted Rt- 
The torque associated with disk truncation and infall of 
material from Rt to the stellar surface is the accretion 
torque Tq. There is also a magnetic torque associated 
with the magnetic connection to a range of radii in the 
disk. This magnetic torque is denoted Tm- The method 
and equations used for calculating i?t, Ta, and Tm are 
given in Paper I and are not modified for the present 
work. 

Two key parameters involved in the calculation of the 
star-disk interaction torques capture the physics of the 
magnetic coupling strength (parameter /3) and the open- 
ing of magnetic field lines due to the differential ro- 
tation (parameter 7c; see Paper I and [Matt fc Pudrit^ 
l2005bf) . Parameter 7c reprents the maximum ratio be- 
tween the azimuthal and vertical magnetic field compo- 
nents {B^/Bz) threading the disk, in order for the dipo- 
lar magnetic field lines to remain closed. In the present 
work, we a dopt 7c = 1 which cor responds to realistic ex- 
pecations (Uzdens kv et al.l l2002^. Parameter 13 is equiv- 
alent to the inverse of the magnetic Reynolds number in 
the inner disk re gion. Here we adopt /3 = 0.01, presented 
in Paper I (and lMatt fc Pudritd l2005b| ) as a reasonable 
guess for T Tauri systems. The effect of varying these 
two parameters was the subject of Paper I. 

2.4.2. Stellar Wind Torque 

In order to calculate the angular momentum loss via 
stellar winds, we use the stellar wind torque formulation 
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Figure 1. Time-evolution of various quantities for the B* = case. Panel (a) shows the rotation period in days, and panel (b) shows 
the spin rate expressed as a fraction of breakup speed. The four lines shown represent the two choices of accretion history and two initial 
spin rates. The red and black lines correspond to the high and low accretion rates, respectively. The solid and dashed lines correspond 
to m odels with an initial spin rate of 0.06 and 0.3 times the breakup rate, respectively. Panel (c) shows the net torque on the star (eq. 
1101 ) for each case. In this panel, the dashed and solid lines of a given color lie on top of each other, since the spin rate does not affect the 
torque, when = 0. The horizontal dotted line indicates T* = for clarity. Panel (d) shows the location of the truncation radius {Rt , 
thick lines) and the corotation radius {Rco, thin lines; eq. [7]), in units of the stellar radius. Here, Rt = -R* for all models, since there is no 
magnetic field to truncate the disk above the stellar surface, so all thick lines (solid and dashed, black and red) overlap. 



of iMatt fc Pudritd (|2008af l. This form ulation is based 
upon analytic calculations going back to iWeber fc Davis! 
([T967), in which the torque from a one-dimensional wind 
is given by 

T„ = ~M.^n,rl, (8) 

where ta is the Alfven radius, the location where the 
wind speed equals that of magnetic Alfven waves. In 
real stellar winds, the flow is as pherical, due to ( e .g.) ro- 
tation and magnetic geometry. IMatt fc Pudritd (|2008al) 
computed solutions for two-dimensional (axisymmetric) 
solar-like stellar winds, using numerical magnetohydro- 
dynamic simulations. They used the simulations to com- 
pute T„ and and then used equation ([5]) to define ta, 
which in this context represents the mass-loss-weighted 
average of the Alfven radius in the multi-dimensional 
flow. For variations in the magnetic field strength, stellar 
radius, surface gravity, and mass loss rate, the solutions 
for TA are well-fit by the following 



(|2008a|) indicate K ss 2.11 and m w 0.223, which we 
adopt in this work. A similar numeri cal study of the 
torque from the winds of massive stars ()Ud-Doula et al.] 
120081 l2009f ) found a relationship similar to equation (jO]) , 
with m = 1/4 and K of the order of unity. By contrast, 
studies of cold, magnetocentrifugally-d r iven flows (e.g., 
Pelletier fc PudritdTl99l iSpruitI 119961 : lAnderson etall 



2005D imply m = 1/3. 



2.4.3. Total Torque and EqmUbrium Spin Rate 

The total torque on the star, used in equation ([5]), is the 
sum of the star-disk interaction torques and the stellar 
wind torque, 



(10) 



The spin rate of any star-disk interacting system gener- 
ally evolves in such a way as to tend toward a theoretical 
"equilibrium spin rate," at which the net torque on the 
star is zero: 



mm 



T* = 



(in spin equilibrium). 



(11) 



(9) 



where K and m are fit constants. For the ca ses with a 
dipole magnetic field, the simulations of lMatt fc Pudrita 



For the values of the magnetic field opening and cou- 
pling par ameters considered in th is paper (7c = 1 and 
/3 = 0.01). IMatt fc Pudritzl ()2008bD showed that the mag- 
netic torque Tm is generally negligible compared to the 
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Figure 2. Results for B» = 500 G and Mm/ Ma = 0.1, in t he s ame format as Figure [T] In panels (a) and (b), the dotted red and 
black lines show the theoretical equilibrium rotation rates (eq. 1121 1 for the high and low accretion rates, respectively. In panel (d), the 
dash-triple-dotted red and black lines show the stellar wind Alfven radius ( ^2.4.21 for the high and low accretion rates, respectively. Also 
in panel (d), the thin and thick solid lines of a given color and type are indistinguishable, since Rt is very close to -Rco in all four cases. 

stellar wind torque T^. Thus, the spin equilibrium is 
approximately characterized by Tq = — T^,. This pre- 
dicts an equilibrium spin rate, express ed as a fraction 
of bre akup speed, given by (eq. [18] of iMatt fc Pudritzl 
l2008bD . 



Table 1 

Model Parameters 



cq,sw 



= ii:-3/2^(6>n-3)/4^-3m/2^ 



where 



^csc 



(12) 
(13) 
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Equation (jl2p assumes that the truncation radius is al- 
ways near the corotation radius, which turns out to be 
true in the models with non-zero magnetic field presented 
in this paper. Below, we will use the predicted equi- 
librium spin rate of equation (|12|) to compare with the 
results of our spin evolution models. 

2.5. Numerical Method 

The coupled equations ([3]), ([4]), and ([5]) describe the 
evolution of the system in time. We wrote a computa- 
tional code that solves these simult aneously, us i ng th e 
fourth-order Runge-Kutta scheme of iPress et all ()1994[ ). 
starting from = 3 x 10^ yr and ending at 3 Myr. The 
code is described in more detail in Paper 1. 

3. RESULTS 



Table[l]contains the parameters for each case presented 
in this section, listed in order of their presentation and 
grouped by the figures in which the results appear. For 
each case, we present the evolution of the system for the 4 
possible combinations of 2 different initial mass accretion 
rates (parameterized by Mao) and 2 different initial spin 
rates (/o). The "i?* = 0" case is identical to that of 
Paper 1, and we include it here to facilitate comparison 
with the new results (cases W1-W3). 

It is instructive first to examine the non-magnetic case, 
= 0, in which there are no magnetic fields. Figure [T] 
illustrates the evolution of the spin rate — expressed as 
the spin period (panel (a)) and as a fraction of breakup 
speed (panel (b)) — the total torque experienced by the 
star (panel (c)), and the radial locations of the disk inner 
edge [Rt) and corotation radius (i?co) (panel (d)). In the 
= case, the disk extends all the way to the surface 
of the star [Rt = i?*, as indicated by the thick lines in 
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Figure 3. Results for _B» = 2000 G and M^j/Ma = 0.1, in the same format as Figure[2] In panel (d), the thin and thick lines of a given 
color and type are indistinguishable, since Rt is very close to Rco in all four cases. 



panel d), and the torque experienced by the star (panel c) 
is solely the accretion torque, Ta — Ma^/GM^Rt (since 

'^w 0)- 

In all 4 models shown in Figure [U the stellar spin rate 
increases monotonically in time, due to the accretion of 
disk material with high specific angular momentum and 
due to the contraction of the star. Only the most extreme 
model (with the slowest initial spin and lowest accretion 
rate, black solid line in panels (a) and (b)) has a spin 
period in the range 1-3 days during the age range of 1-3 
Myr. The other 3 models have spin periods substantially 
less than a day, in the same age range. 

All subsequent cases (W1-W3) have a non-zero mag- 
netic field and use 7c = 1 and (3 — 0.01 for the star-disk 
magnetic connection parameters, which are appropriate 
values to include the effects of strong magnetic coupling 
to the disk and the op ening of magnetic field lines by 
azimuthal twisting fsee lMatt fc Pudrit3l2005bl and Pa- 
per 1). Thus, in addition to the accretion torque, the 
model stars experience spin-up and spin-down torques 
(Tm) associated with a magnetic connection over a fi- 
nite region of the disk. Also, cases W1-W3 include the 
mass loss (x) and torque (T„,) due to stellar winds (as 
described in Sj2]). 

Figure [2] shows the results for the four models in the 
Wl case, in which the star has a 500 Gauss magnetic 
field and stellar wind mass loss rate that is 10% of the 
disk accretion rate (x = 0.1). In this Figure, the line 
colors and styles have the same meaning as Figure [U 



with a few additions. Panels (a) and (b) include dotted 
lines that indicate the equilibrium spin value predicted 
by equation (fT2|) for the high and low mass accretion rate 
models (red and black lines, respectively). In panel (d), 
the dash-triple-dotted lines indicate the location of the 
Alfven radius in the stellar wind (eq. [9] ) for the models 
with high and low accretion rate (and thus high and low 
wind outflow rate; red and black lines, respectively). 

By comparing Figures [T] and [21 it is clear that the stel- 
lar spin evolution is substantially influenced by the pres- 
ence of a 500 G held and accretion-powered stellar wind 
with X = 0.1. In the age range of 1-3 Myr, all four mod- 
els in Figure[2]have spin periods in the range of 1-5 days, 
which corresponds to spin rates between 10% and 20% 
of breakup. As visible in panel (d), the truncation ra- 
dius is very close to the corotation radius (such that the 
thick and thin lines overlap in the plot), for all models 
at all times. The Alfven radius in the stellar wind (panel 
(d)) for all models remains in the range of approximately 
10-30 i?*. 

Panels (a) and (b) show the equilibrium spin rate (dot- 
ted lines) predicted by equation ()12p under the assump- 
tion that the spin evolution is dominated by a balance 
between the accretion (spin-up) torque and stellar wind 
(spin-down) torque. In Figure [2j the two models corre- 
sponding to a high accretion rate (red solid and dashed 
lines) approach the equilibrium spin rate at an age of 
a few times 10^ yr. After this time, the initial condi- 
tion for the spin rate of these models has effectively been 
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Figure 4. Results for B, = 2000 G and M^j/Ma = 0.01, in the same format as Figure[2l In panel (d), the thin and thick lines of a given 
color and type are indistinguishable, since Rt is very close to Rco in all four cases. 



"erased," in a sense that the subsequent spin rate is in- 
sensitive to the initial rate. For the models with low 
accretion rate (black solid and dashed lines), the torque 
is not strong enough to drive the spin rate to the equi- 
librium value within 3 Myr. The case with a slow initial 
spin rate (black solid line) evolves near the equilibrium 
value, but this is only due to a coincidence between the 
initial and equilibrium spin rates. 

Figure |3] shows results for the W2 case, which is the 
same as Wl, except that the magnetic field is 4 times 
stronger (B* = 2000 G). By comparing Figures [2] and 
[3l it is clear that the stronger magnetic field results in 
slower spin rates. In the age range of 1-3 Myr, all four 
models in Figure |3] have spin periods in the range of 4- 
12 days, corresponding to spin rates between 3% and 8% 
of breakup. The models with high accretion rate in the 
W2 case (Fig. |31 red solid and dashed lines) approach 
the equilibrium spin rate after ^ 10^ yr, sooner than 
the weaker field case (Wl). Furthermore, the models 
with low accretion rate (black solid and dashed lines) 
approach the equilibrium spin rate after ^ 10^ yr. In 
all models and at all times, the truncation radius is very 
close to the corotation radius (such that the thick and 
thin lines overlap in panel (d)), and the Alfven radius in 
the winds remains in the range of 15-50 i?*. 

Figure S] shows results for the W3 case, which is the 
same as W2, except that the stellar wind outflow rate 
is 10 times less (x — 0.01). A lower wind mass outflow 
rate, if all else is equal, means a lower spin-down torque 



from the stellar wind. Thus, the accretion-powered stel- 
lar wind in the W3 case is less effective at spinning down 
the star than the W2 case. A comparison between Fig- 
ures |4] and [2] reveals that the W3 case is qualitatively and 
quantitatively very similar to the Wl case, which has a 
higher wind outflow rate but a weaker magnetic fleld. 
The similarity in the spin evolution of these two cases is a 
consequence of the dependence of the stellar wind torque 
on the parameters changed. Specifically, equations ([5]) 
and © and TO = 0.223 resufls in T„ cx M°-^^B°-^^, 
which means that the factor of 10 difference in is 
almost completely compensated by the factor of 4 differ- 
ence in B^. While the stellar wind torques are similar 
in these two cases, the Alfven radius in the wind has a 
different dependence on these parameters, and thus the 
models in the W3 case have a much larger Alfven radii 
than the Wl case (compare the dash-triple-dotted lines 
in panel (d)). 

It is instructive to compare the cases presented here 
with the cases in Paper 1. Case Wl is identical to case 
01 in Paper 1, except that Wl includes the effect of an 
accretion-powered stellar wind, via a nonzero value of the 
mass-loss parameter x- Likewise, cases W2 and W3 are 
similar to 02 of Paper I. A comparison between W1-W3 
(Figs, in and 01/02 (Figs. 6 and 7 of Paper I) indi- 
cates that the cases with stellar winds spin substantially 
slower than the cases without. For the range of magnetic 
field strengths, accretion rates, initial spin rates, and 
mass outflow rates considered in cases W1-W3, the mod- 
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els exhibit rotation periods within the range of approx- 
imately 1-10 days in the age range of 1-3 Myr. These 
rotation periods correspond to spin rates in the range of 
approximately 5%-20% of breakup speed and lie within 
the bulk of the distribution observed in stars of similar 
age and mass (see SJT|). 

For all of the models in cases W1-W3 that approach an 
equilibrium, the spin rates closely follow the prediction 
of equation (IT2]) . Since that equation assumes a balance 
solely between the accretion torque and a stellar wind 
torque (T, = —T^), the fact that some models follow this 
prediction indicates that the magnetic torque from the 
star-disk conne ction (Tm) is negligible in those models. 
More generally. iMatt fc Pudr"iti(|2008b[ ) showed that, for 
the strong coupling case considered here (i.e., (3 = 0.01 
and 7c = 1), the stellar wind torque dominates over the 
spin-down portion of the magnetic torque Tm as long as 
the Alfven radius in the stellar wind (ta) is < 84i?*. 
This is true during the entire evolution for all models 
presented here (see dash-triple-dotted lines in panel (d) 
of Figs. [2H4]). Thus, it is the stellar wind torque, not the 
star-disk magnetic connection, that is responsible for the 
slower spin of these models relative to the non-magnetic 
case. 



1.00 



4. COMPARISON BETWEEN APSW AND DISK LOCKING 
MODELS 

For the magnetic torque model and chosen parameters 
presented in sections [2] and |3l we showed that the spin- 
down torque from APSWs completely dominates over 
any spin-down torques that may arise from the mag- 
netic connection between the star and disk. Thus, the 
evolution of the stellar spin is primarily determined by 
a competition between the spin-down torques from an 
APSW and a spin- up via the accretion torque (hereafter, 
the "APSW model"). However, given the large number 
of unknowns regarding real systems, it is instructive to 
compare the predictions of different star-disk interaction 
models for a range of possible conditions and assump- 
tions. 

As discussed earlier, in these systems the spin rate of 
the star always tends toward the equilibrium spin rate at 
which the net torque on the star would be zero. If the 
spin- up or spin-down time is short (i.e., for a strong net 
torque) compared to the evolution timescale, the star's 
spin will remain close to this equilibrium value. Whether 
or not the torque is strong enough to keep the star spin- 
ning near equilibrium depends upon many factors (e.g., 
the accretion rate Ma and magnetic field strength i?,), 
which also may be varying in time. Thus it is not gener- 
ally expected that all systems should be in equilibrium. 
However, for the purpose of comparing different mod- 
els, it is convenient and particularly instructive to com- 
pare the predictions of the equilibrium spin rate for each 
model. When a particular model predicts a slower spin 
rate compared to other models, it simply means that the 
spin-down torques are generally stronger in that model. 

The equilibrium spin rate predicted by the APSW 
model, expressed as a fraction of breakup speed, is given 
in equation ((T2|) . In Figure [3 the solid lines show this 
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Figure 5. Predicted equilibrium spin rates of different star-disk 
interact ion models as a function of the dimensionless parameter 
'3/ (eq. I13I V The solid lin es show the prediction of the APSW 
model (determined by eq. 1121 . with K = 2.11 and m = 0.223) 
for three different values of Xi indicated. The dashed lines 
show the predictions for t he Gh osh &; Lamb- type dis k l ocking 
model of IMatt fc Pudritj I l2005bl ') (determined by eq. [M], with 
C = C{l3,jc)) for 7c = 1 and two different values of the mag- 
netic coupling parameter /3, as indicated. The dotted li ne s hows 
the prediction of the X-wind model (determined by eq. 1141 with 
C = 0.974). 



predicted equilibrium spin rate as a function of the di- 
mensionless parameter 4' and for 3 different values of 
Mw/Ma = X — 1/3, 0.1, and 0.01. As discussed in sec- 
tion [21 we adopt values oi K = 2.11 and m = 0.223, 
which results in a power-law slope of — 3m/2 —0.33 
in the Figure. The plot covers a range in "i! from 1 to 
10^, which is chosen to approximately reflect the possible 
range in real systems, obtained by choosing extreme val- 
ues of parameter s from the range that is observed for T 
Taur i stars (e.g., iHartmann et al. 119981: iGullbring et al.l 



fl^.'Johns-KruU fc Gafford"20 02t iRobberto et aLl l2004l: 
Siciha-Aguilar et al. 2006; Bouvi er et al.ll2007[ ). 

Of the various star-disk interaction studies in the 
literature, there are two other types of general mod- 
els that predict an equilibrium spin rate, and with 
which we will compare the predictions of the APSW 
model. These are the Gh osh &: Lamb type (GL- 
type) models (e.g., Paper I; 'Camenzind"1990i 'Koenigl 
[1991; Cameron & CamobeU 1993; Cameron et al. 1995[ 
Ti 1994, 1995; Lovelace et al. 1995; Armitage & ClarS 
,1996; Matt fc Pudrita.2005b : .Zweibel et al...2006 ). which 
are base d on the original ac creting magnetized star 
model of|G hosh fc Lam bI (119781). and the X-wind model 
(IShu et aLr:1 994a.b: N aiita fc Shulll99l lOstriker fc Shul 
119951 : IShu et al.lll99,5t IMohantv fc Shdl2008t) . Both the 
GL-type and X-wind models assume that the magnetic 
connection between the star and disk alone is responsible 
for removing the necessary angular momentum, and both 
neglect the influence of stellar winds. In their equilibrium 
state, t hese are often referred to as "disk locking" models 
(|Choi fc Herbst,,1996i) . The predicted equilibrium spin 
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rate for both types of models can be expressed as 

/cq41 = C'f-''/\ (14) 

where C is a dimensionless constant that takes into ac- 
count various assumptions in the models, and ^' is the 
same as defined in equation (jl3p . 

The various disk-locking models in the literature dif- 
fer in their assumptions, which is reflected in a dif- 
ferent value of C. T he star-disk torque model of 
IMatt fc Pudritd (|2Q05bf ) (used in the present paper) is 
a GL-type model that includes the effect of magnetic 
field opening (via the parameter 7c) and magnetic cou- 
pling strength (via the parameter /3). That work showed 
that the constaniQ C is a function of only jc and /?, and 
different choices of these can result in C values that span 
the range of GL - type models presented in the literature. 
lUzdenskv et al.l (|2002l ) demonstrated that a value 7c = 1 
(adoped in the present work) is the most appropriate to 
properly take into account the effect of the opening of 
magnetic field lines due to differential twisting between 
the star and disk. 

The dashed lines in Figure [5] show the prediction of 
GL-type models for 7c = 1 and two different values of 
P (0.01 and 1.0, corresponding to C ~ 10.1 and C w 
2.49, respectively). The upper-most line corresponds to 
the GL-type model used to compute the magnetic star- 
disk torques in the spin-calculations of the present paper 
(i.e., for P = 0.01). It is clear that this GL-type model 
predicts an equilibrium spin rate that is faster than the 
APSW model, for the entire plotted range of ^, and as 
long as X is greater than a few thousanths. This is just 
another way of demonstrating that, when the mass loss is 
high enough, APSWs are much more effective at spinning 
down stars than a GL-type star-disk interaction alone. 

The effective value of /3 in real systems is not well 
constrained. In the G L-type models with 7c = 1, 
IMatt fc Pudritd (j2005b( ) showed that the value of /3 = 1 
results in the minimum value for C and thus represents 
the case where the GL-type model torques would be the 
most effective at spinning the star down. The bottom- 
most dashed line in Figure [5] shows the prediction for this 
special case. It is clear that the APSW model generally 
predicts a slower spin rate for x ~ 0.1 than the GL-type 
models with any valu£0 of /?. For x 0.01, the APSW 
model predictions are of a similar magnitude as those of 
a GL-type model with (3 — 0.1. While the effective value 
of (3 in real systems is uncertain, it seems likely that the 
inner disks of young stars exist in a state of high magnetic 
Reynolds number {/3 ^ 1), and thus that the GL-type 
spin-down torques are relatively weak. 

In the X-wind model, the angular momenum is as- 
sumed not to accrete onto the star, but rather is inter- 
cepted at the inner edge of the disk (the "X-point" ) and 
then carried out of the system by a wind (the "X-wind" ) 
from the X-point. The fiducial model of iShu et all 

^ For convenience here, we have absorbed all constant factors 
into t he value of C in e quation II14I I. as compared to equation (25) 
of , Matt fc Pudritd pOOSbV 

^ This statement is true for 7c = 1. Most of the GL-type mod- 
els in the literature have not taken into account the opening of 
the magnetic field, which is equivalent to assuming 7c = oo- In 
this case, an d especially when /3 is take n to be small (e.g., as in 
the model of lCameron iz CamDbet]||1993l ). the spin-down torque is 
(artificially) enhanced and can be arbitrarily large. 



(|1994aD and lOstriker fc Shul ([l995l) supposes that the 
wind outflow rate from the X-point is ~ 1/3 times the 
accretion rate, which results in a predicted equilibrium 
spin rate that is given by equation ()14|) . with a value of 
C 0.974. The dotted line in Figure [5] shows this pre- 
dicted spin rate of the X-wind model. It is clear that this 
model also predicts slower spin rates than the GL-type 
models. 

Both the APSW and X-wind models employ a wind 
to remove angular momentum from the system. Thus, 
it is perhaps not surprising that the predicted spin rate 
of the APSW model with the same mass outflow rate 
as the X-wind, x = 1/3 (bottom solid line in Fig. [S]) , is 
quantitatively similar to the predictions of the X-wind. 
Furthermore, it can be shown that the factor C in the 
X-wind model depends upon the mass loss rate to the 
power of —3/7, which is very close to the dependence of 
the APSW model prediction on the mass loss rate (in eq. 
[12], / cx x^O '*^, for the adopted value of m « 0.223). 

Although the predicted equilibrium spin rates are 
quantitatively similar, the physical picture of APSW and 
X-wind models have some important differences. In the 
X-wind, the wind is magnetically connected to, and re- 
moves angular momentum from, the X-point (not from 
the star, as in a stellar wind). The X-wind model is de- 
veloped under the assumption that there is no net torque 
on the star and that the X-point is located at the corota- 
tion radius. Thus, the balancing of angular momentum 
is evaluated at the X-point (rather than on the star), and 
this leads to a requirement for a certain amount of stellar 
magnetic flux to be pinched, or "trapped," within the X- 
point. Having the flux is trapped at the X-point avoids 
the large-scale opening of magnetic field (which is a prob- 
lem for the GL-type models) and also produces a ge- 
ometry in which magneto-centrifugal forces may provide 
a natural explanation for the acceleration of the wind. 
However, the biggest problem for the X-wind model is 
that the amount of flux trapping required at t he X-point 
raise s questions of dynamical stability (e.g., lUzdenskyI 
12004 ), and such a magnetic field configuration has not yet 
been demonstrated to exist in a dynamical model. Fur- 
thermore, since the model assumes no net torque on the 
star (i.e., it assumes the star to be in spin equilibrium), 
the model provides little information about how angular 
momentum may fiow from accreting material outward to 
the X-point, and accross the X-point, in order to be re- 
moved in the X-wind and/or via transport processes in 
the disk ()Ferreira et al.ll2000l ). The assumption of equi- 
librium also means that the model cannot address how 
these systems evolve (e.g., one cannot make plots such as 
panels [c] in Figs.[TH4]), nor whether or how a system will 
reach or approximate the equilibrium state and magnetic 
field configuration described. 

By contrast, the APSW model calculates the spin- up 
and spin-down torques acting on the star from various 
processes (i.e., at different geographical locations on the 
star). The theoretical equilibrium state corresponds to 
when the net torque is zero, but the model can be used 
to compute the net torque in any state. The magnetic 
field configuration is that of the stellar magnetic field, 
which has a significant amount of fiux that is opened 
by various dynami c al pro cesses (e.g., see discussion 
in IMatt fc Pudritd l2005b( ). Numerical, dynamical 
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20 lit) generally show that such a magnetic field con- 
figuration is a natural and general consequence of the 
star-disk interaction. It is the open stellar field lines 
that are primarily responsible for the spin-down torque 
on the star, via a wind flowing along the field. The 
amount of magnetic flux participating in the wind is not 
a free parameter, but is determined in a dynamically 
self-consistent way via the MHD simulations used 
to compute equation (O for the stellar wind torque 
(jMatt k Pudritd 12008a). Thus, there is no conceptual 
nor dynamical problem with the magnetic field configu- 
ration of the APSW model. However, the biggest open 
question for APSW model is whether and/or how the 
accretion power in the system can drive a wind that is 
magnetically connected to the star and that has a mass 
outflow rate high enough to extract significant angular 
momentum. This is discussed further in section O 

It is clear from Figure [5] that the APSW predicts 
a slope that is somewhat shallower than the disk- 
locking models. This fact could in principle be used 
to discriminate between th e models (e.g., similar to 
iJohns-Kruh fc GaffordI [20021. ahhough this may not be 
possible with existing data. The difficulty is due to at 
least 3 factors: a) there are large observational uncer- 
tainties for the quantities plotted in Figure [S] and few 
stars for which these quantities have been measured; b) 
not all systems are expected to be in spin-equilibrium; 
and c) these systems are typically highly time-variable 
(e.g.. lHartmanniri997l ). which means that a realistic spin- 
equilibrium (net zero torque) state may only exist in a 
time-averaged sense — or in a statistical sense, for a large 
collection of stars. Furthermore, the predicted slope for 
the APSW model may be different, depending upon a 
few unknowns. The form of equation ()12|) assumes that 
the d isk truncation radius is close to the corotation ra- 
dius ()Matt fc Pudritdr2008b[) . which will not necessarily 
be true in all systems. Also, the value of the power law 
parameter m depends upon the conditions in the stel- 
lar wind, and a number of factors (discussed in Sj5]) may 
affect this value. 

Finally, we note that in the interest of performing di- 
rect and quantitative comparisons, this section has fo- 
cused on the three classes of star-disk interaction models 
for which a simple prediction for a hypothetical equilib- 
rium spin rate is possible. Thus, there are a number of 
other ideas in the literature that we have neglected and 
that nonetheless deserve further attention. For example, 
time- variability of the accretion rate onto the star may 
stron gly influence stellar spiii evolution (e.g., iPophaml 
119961 : iD'Angelo fc Sprui"tll2010D . and the dynamics of the 
interface region that exists between the open magnetic 
flux threading the star and disk may also be i mpor- 
tant for extracting angular momentum (e.g ., Hiros e et al.l 
Il997t iFerreira et all 120001: iRomanova et al.M2009. Zanni 
& Ferreira, in preparation). 

5. SUMMARY AND DISCUSSION 



Paper I presented a model for computing the evolution 
of mass, radius, and spin rate of an accreting, one solar 
mass star during the Hayashi phase (from 3 x 10** yr 
to 3 Myr). The spin of the star evolves from its initial 
value due to changes in the moment of inertia (due to 
stellar contraction) and to external torques. To compute 
the external torques. Paper I included only the magnetic 
star-disk interaction torque formulation of MP05, and 
neglected any torque along open field lines. That work 
showed that, for strong magnetic coupling to the disk 
(7c = 1 and (3 ~ O.Ol), cases with either i?, ~ 500G 
or 2000G result in spin evolution that is only slightly 
changed from the non-magnetic case. In other words, 
when the magnetic coupling is strong and -B* < 2000G, 
the spin-down torque arising in the star-disk interaction 
is negligible for the spin rates and mass accretion rates 
considered. 

The goal of the present paper was to determine under 
what conditions the additional spin-down torque from 
stellar winds may significantly affect the spin evolution 
and produce a range of spin rates consistent with the 
observed range. Thus, we have extended the model of 
Paper I (as described in to include an additional 
spin-down torque due to the stellar wind. As stellar 
wind torques depend upon most of the same parame- 
ters as star-disk interaction torques, including a stellar 
wind has introduced only one new parameter, the mass 
outflow rate of the wind. Here, we consider the case 
where the mass outflow rate is some fixed fraction (x) of 
the accretion rate, as may be expected for a wind that 
is som ehow driven by accretion power ()Matt fc Pudritj 
l2005aD . 

In section [3l we found that the cases with B = 2000G 
and mass loss rates of x = 0.1 resulted in spin rates 
in the range of ~ 5-10 days, in the age range of 1-3 
Myr. For the two cases with either a lower mass loss rate 
(x = 0.01) or magnetic field strength {B = 500G), the 
resulting spin rates were in the range of 1-5 days (for 
1-3 Myr ages). These two cases exhibited similar spin 
evolution because the stellar wind torques were nearly 
identical in both cases (i.e., there is a degeneracy be- 
tween and in the stellar wind torque, eq. [8] and 
[S]; discussed in Overall, the magnetic cases pre- 

sented in section |3] resulted in stellar spin periods in the 
range of approximately 1-10 days, in the age range of 
1-3 Myr. This range covers the range of the most pop- 
ulated regions of the observed spin period dis tributions 
of clusters with similar ages (e.g., lStassun~et al. 19991 
iHerb st et al. | l2002l: [Rebuh et al.ll2004l: iLamm et al.ii200£ 



Herb st et all 120071: Ih^in fc Bouvie'Jl2009[) . In general, 
accretion-powered stellar winds may explain the slowest 
rotators as those stars with the strongest magnetic fields, 
lowest accretion rates, and/or highest stellar wind mass 
outflow rates (relative to the accretion rate) . 

We found that some cases that had torques strong 
enough to drive the system close to an equilibrium spin 
rate (e.g., as in Fig. [3l). Spin equilibrium states are in- 
teresting because the initial conditions are "erased," and 
the spin rate depends only on present conditions (e.g., ac- 
cretion rate, magnetic fleld strength, etc.). However, for 
the conditions considered here (in particular, for strong 
magnetic coupling), the magnetic spin-down torques aris- 
ing in the star-disk interaction were completely negligi- 
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ble. Thus, these stars' equilibrium rotation rates are not 
physically linked to the rotation rate of the disk inner 
edge (as in disk locking) but rather are described by a 
simple balance between a spin-up torque from accretion 
and a spin-down torque from a stellar wind. This simple 
balance allows for an analytic formula for predicting the 
equilibrium (net-zero torque) spin rate (e.g., eq. [H]). 

While we have necessarily made a number of simpli- 
fying assumptions, we have whenever possible adopted 
the same assumptions that are typically used in previous 
models in the literature (in particular, the disk-locking 
models). Our approach in Paper I and here has been 
deliberately systematic, so that each additional compo- 
nent or assumption in the model can be understood in 
turn. In this way, our approach best serves to highlight 
the influence of the opening of the mangetic field due 
to star-disk differential rotation (Paper I) and the influ- 
ence of an acctetion-powered stellar wind (present pa- 
per), relative to the previous studies that do not include 
these effects. We have explored a range of conditions 
that are consistent with observations of T Tauri stars, 
but we have not explored all of parameter space, and 
there remain uncertainties inherant to various adopted 
assumptions and approximations. Thus, while our quan- 
titative results should be viewed as approximate, the re- 
sults relative to other studies/models in the literature are 
robust. 

With this in mind, we compared the predicted equilib- 
rium spin rate of the APSW scenario to the predictions 
of two types of disk-locking models, the Ghosh & Lamb 
type and X-wind models O. Overall, this comparison 
and our results in general demonstrate that APSWs can 
explain the observed distribution of young star spins in 
a similar way as the classical disk locking picture, while 
at the same time avoiding the problem of magnetic fleld 
line opening (for the GL-type models), as well as the 
assumption of spin equilibrium and requirement of sig- 
nificant flux trapping (for the X-wind model). 

The APSW scenario has one additional parameter, the 
mass loss rate in the wind. For the assumptions adopted 
here, we found that in order for the APSWs to have 
a significant infiuence, the mass loss rates should be at 
least of the order of a percent of the accretion rate. It 
still remains to be shown whether and/or how the en- 
ergy derived from the accretion process may drive a wind 
that is magnetically connected to the star, and with suf- 
ficiently high mass outflow rate. Such winds would not 
likely be driven signiflc antly by thermal press ure (due 
to a rapid cooling time iMatt fc Pudritd I2007D , nor by 
magneto-centrifugal effects (from slowly -rotating stars) , 
but Alfven waves may be important (e.g.,'Decampli"1981l; 
Hartmann et al. 1982; Suzuki 2007). Cranmcr (2008, 
2009( 1 demonstrated that Alfven waves generated by the 
accretion process in T Tauri stars are capable of driv- 
ing enhanced stellar winds. The mass loss rates de- 
rived by those models were typically x ^ 0.01. There 
is a hard energetic upper li mit of % < 0.6 for AP- 
SWs (IMatt fc Pudritzl l2008b[ ). and values close to this 
will have observational conseq uences that may a.lread y 
be ruled out in some systems (jZanni fc Ferreirall2011[ ). 
Thus, values of x much greater than 10-20% seem un- 
likely, in general. Clearly, more work is needed to deter- 
mine what is the mass loss rate along the stellar magnetic 
field and how this should depend upon system parame- 



ters. 

Measurements of the magnetic field strengths of young 
stars suggest that the stellar surface is blanketed with 
comp lex magnet i c fields with a strengths of ^ 2 kG 
(e.g,.iS afiCT 1998; Bouvi er et all [20(171 : 1.Tohns-KrullI [20071 : 
lYang fc Johns-KruU 20^^ However, even when com- 
plex magnetic fields are present, it is the large-scale 
(dipole) component that is the most important for 
the star-disk interac tion and torques on the star (e.g., 
iGregorv et al.l 120081 ). The dipole components have 
been measured via spectropolarim etry for at le ast 10 
accre ting T Tauri stars to date (jP onati et al.' 200% 
2008tlHussain et al.ll2009 HDonati et al...2qi0a.b, . ,2011bjal : 
Donati. J. F. et al. 2011: iSkellv. M. et al.ll2(ni[ ). and the 
equatorial field strengths of the dipole components range 
from ~ 100 — 1000 G, with the top of this range deter- 
mined by just two of the stars (i?* ~ 600 G for BP Tau , 
and w 1000-1500 G for AA Tau ; iDonati et al1[200l 
l2010aD . The disk- locking models usually require field 
strengths of ^ lOOOG to explain the slowest rotators, and 
the weakest equatorial field strength we have considered 
in the present work is B^, = 500G. The APSW scenario 
can, in principle, make up for a weaker magnetic field by 
having a larger mass outfiow rate, although the largest 
values of x considered in the present p aper are already 
approaching the upper limits (see, e.g.. IMatt fc Pudritd 
I2008bl: [Zanni fc Ferreirall201lD . To better constrain the 
models, it will be important to have magnetic field mea- 
surements, particularly of the large-scale components, for 
a larger sample of pre-main-sequence stars. 

In addition to those discussed above, there remain a 
number of caveats to the present work. While we have 
adopted a relatively sophisticated model for the torques 
on the star, we use a simplified treatment for the accre- 
tion history, the evolution of the magnetic field, and the 
structure and evolution of the star itself. In order to be 
able to make a more meaningful comparison with obser- 
vations, and especially to be able to evolve the system 
to much later times, it will be necessary to improve the 
theoretical treatment of these components. 

For example, the present treatment of the accre- 
tion does not allow for a possible state of the system 
in which the disk is truncated outside of the corota- 
tion radius, the so called "propeller" regime. In this 
regime, there is generally no accretion onto the star (e.g . , 
llUarionov fc SunvaevI 119751 : iSunvaev fc Sh akuTa T977|, 
or the accretion is intermitte nt (e.g., iRpmanova et aj] 
l2005HD'Angelo fc Sprui"tll2010( ). In the present work, we 
are considering systems with strong magnetic coupling 
(large magnetic Prandtl number in the disk) and rela- 
tively slow rotation. Under these conditions, the mag- 
netic spin-down torque on the star — which acts to spin 
up the disk and can potentially truncate the disk outside 
of Rco — is relatively weak compared to other torques in 
the system. Thus, we don't expect the propeller regime 
to be important for most of the evolution these systems. 
However, there may be times in the history of such stars 
(e.g., toward the end of the accretion phase) in which 
the propeller regime is important for the angular mo- 
mentum loss. Properly capturing the transition between 
accreting and propeller states requires a self-consistent 
trea tment of the evolution of the accretion disk (as in , 
e.g.. lArmitage fc Clarkel[T996t ID'Angelo fc Spruitl[20Toh . 
In the present model, we imposed the accretion rate onto 
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the star, which implicitly assumes that the disk is always 
able to transport ex cess angular momentum given to it 
by the star (see, e.g. JMatt fc Pudritd[2005bl ). To explore 
the influence of the propeller regime on the long-term 
spin evolution, future models should include both self- 
consistent disk evolution and sophisticated treatment of 
the torques on the star. 

Finally, all of our conclusions regarding the effective- 
ness of stellar wind torques re hes on the particular for - 
mulation for the torque from IMatt fc Pudrit j (|2008aD . 
While this is the most appropriate formulation for low- 
mass-star winds that exists in the literature, there re- 
main a few open questions about its use in the present 
work. First of all, that formulation was derived for an 
isolated star, and we have simply added this torque to 
our model, which also assumes the presence of an accre- 
tion disk and associated interaction torques in a non-self- 
consistent way. It is not yet clear how the presence of 
a disk will influence the stellar wind, and in particular 
how it may influence the dependence of the stellar wind 
torque on p arameters (e.g., the valu e of m). Also, the for- 
mulation of lMatt fc Pudrita ()2008al ) was based upon sim- 
ulations of a star with a singular spin rate (/ = 0.1) and 
a solar-like wind acceleration mechanism. It is encourag- 
ing th at a study of angular momen tum flow from massive 
stars (|Ud-Doula et al.ll2008l[2009l ). which includes a very 
different wind acceleration mechanism and a range of 
spin rates, found a similar power law relationship for the 
torque, but more studies are warranted. Lastly, all stel- 
lar wind torque formulations in the literature are based 
upon simple (e.g., dipolar) magnetic geometries, while 
we know that T Tauri stars possess complex magnetic 
field structures. In future work, it will be important 
to determine whether the presence of an accretion disk, 
more complex magnetic geometries, or different wind ac- 
celeration mechanisms will significanlty enhance or sup- 
pr ess the stellar wind torq ue relative to the formulation 
of IMatt fc Pudritd ()2008al) adopted here. 

Much work remains to develop the necessary theory, to 
improve the precision and number of observational mea- 
surements of relevant system parameters, and ultimately 
to understand the observed distributions and evolution 
of stellar spin rates, but the idea that powerful stellar 
winds may extract significant angular momentum from 
accreting stars remains a promising scenario. 
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